Steady viscous flow past a circular cylinder

Viscous flow past a circular cylinder becomes unstable around Reynolds number Re = 40. With a numerical technique based on Newton's method and made possible by the use of a supercomputer, steady (but unstable) solutions have been calculated up to Re = 400. It is found that the wake continues to grow in length approximately linearly with Re. However, in conflict with available asymptotic predictions, the width starts to increase very rapidly around Re = 300. All numerical calculations have been performed on the CDC CYBER 205 at the CDC Service Center in Arden Hills, Minnesota

A Numerical Methodology for the Painlevé Equations

The six Painlevé transcendents PI – PVI have both applications and analytic properties that make them stand out from most other classes of special functions. Although they have been the subject of extensive theoretical investigations for about a century, they still have a reputation for being numerically challenging. In particular, their extensive pole fields in the complex plane have often been perceived as ‘numerical mine fields’. In the present work, we note that the Painlevé property in fact provides the opportunity for very fast and accurate numerical solutions throughout such fields. When combining a Taylor/Padé-based ODE initial value solver for the pole fields with a boundary value solver for smooth regions, numerical solutions become available across the full complex plane. We focus here on the numerical methodology, and illustrate it for the PI equation. In later studies, we will concentrate on mathematical aspects of both the PI and the higher Painlevé transcendents

State diagram and the phase transition of pp-bosons in a square bi-partite optical lattice

It is shown that, in a reasonable approximation, the quantum state of $p$-bosons in a bi-partite square two-dimensional optical lattice is governed by the nonlinear boson model describing tunneling of \textit{boson pairs} between two orthogonal degenerate quasi momenta on the edge of the first Brillouin zone. The interplay between the lattice anisotropy and the atomic interactions leads to the second-order phase transition between the number-squeezed and coherent phase states of the $p$-bosons. In the isotropic case of the recent experiment, Nature Physicis 7, 147 (2011), the $p$-bosons are in the coherent phase state, where the relative global phase between the two quasi momenta is defined only up to mod($\pi$): $\phi=\pm\pi/2$. The quantum phase diagram of the nonlinear boson model is given.Comment: 15 pages; 5 figures, some in colo

Stable Computations with Flat Radial Basis Functions Using Vector-Valued Rational Approximations

One commonly finds in applications of smooth radial basis functions (RBFs) that scaling the kernels so they are \u27flat\u27 leads to smaller discretization errors. However, the direct numerical approach for computing with flat RBFs (RBF-Direct) is severely ill-conditioned. We present an algorithm for bypassing this ill-conditioning that is based on a new method for rational approximation (RA) of vector-valued analytic functions with the property that all components of the vector share the same singularities. This new algorithm (RBF-RA) is more accurate, robust, and easier to implement than the Contour-Padé method, which is similarly based on vector-valued rational approximation. In contrast to the stable RBF-QR and RBF-GA algorithms, which are based on finding a better conditioned base in the same RBF-space, the new algorithm can be used with any type of smooth radial kernel, and it is also applicable to a wider range of tasks (including calculating Hermite type implicit RBF-FD stencils). We present a series of numerical experiments demonstrating the effectiveness of this new method for computing RBF interpolants in the flat regime. We also demonstrate the flexibility of the method by using it to compute implicit RBF-FD formulas in the flat regime and then using these for solving Poisson\u27s equation in a 3-D spherical shell

Effects of Line-tying on Resistive Tearing Instability in Slab Geometry

The effects of line-tying on resistive tearing instability in slab geometry is studied within the framework of reduced magnetohydrodynamics (RMHD).\citep{KadomtsevP1974,Strauss1976} It is found that line-tying has a stabilizing effect. The tearing mode is stabilized when the system length $L$ is shorter than a critical length $L_{c}$, which is independent of the resistivity $\eta$. When $L$ is not too much longer than $L_{c}$, the growthrate $\gamma$ is proportional to $\eta$ . When $L$ is sufficiently long, the tearing mode scaling $\gamma\sim\eta^{3/5}$ is recovered. The transition from $\gamma\sim\eta$ to $\gamma\sim\eta^{3/5}$ occurs at a transition length $L_{t}\sim\eta^{-2/5}$.Comment: Correct a typ

Spectral methods for the wave equation in second-order form

Current spectral simulations of Einstein's equations require writing the equations in first-order form, potentially introducing instabilities and inefficiencies. We present a new penalty method for pseudo-spectral evolutions of second order in space wave equations. The penalties are constructed as functions of Legendre polynomials and are added to the equations of motion everywhere, not only on the boundaries. Using energy methods, we prove semi-discrete stability of the new method for the scalar wave equation in flat space and show how it can be applied to the scalar wave on a curved background. Numerical results demonstrating stability and convergence for multi-domain second-order scalar wave evolutions are also presented. This work provides a foundation for treating Einstein's equations directly in second-order form by spectral methods.Comment: 16 pages, 5 figure

Anomalous Spin-Related Quantum Phase in Mesoscopic Hole Rings

We have obtained numerically exact results for the spin-related geometric quantum phases that arise in p-type semiconductor ring structures. The interplay between gate-controllable (Rashba) spin splitting and quantum-confinement-induced mixing between hole-spin states causes a much higher sensitivity of magnetoconductance oscillations to external parameters than previously expected. Our results imply a much-enhanced functionality of hole-ring spin-interference devices and shed new light on recent experimental findings.Comment: 6 pages, 4 figures, RevTe